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Abstract 

The problem of finding the asymptotic behavior of the maximal density 0max of sphere packings 
in high Euclidean dimensions is one of the most fascinating and challenging problems in discrete 
geometry. One century ago, Minkowski obtained a rigorous lower bound on 0max that is controlled 
asymptotically by 1/2"^, where d is the Euclidean space dimension. An indication of the difficulty of 
the problem can be garnered from the fact that exponential improvement of Minkowski's bound has 
proved to be elusive, even though existing upper bounds suggest that such improvement should be 
possible. Using a statistical-mechanical procedure to optimize the density associated with a "test" 
pair correlation function and a conjecture concerning the existence of disordered sphere packings 
[S. Torquato and F. H. Stillinger, Experimental Math. 15, 307 (2006)], the putative exponen- 
tial improvement on (/>niax was found with an asymptotic behavior controlled by \^'2^^''^'^^^^"'^^ . 
Using the same methods, we investigate whether this exponential improvement can be further im- 
proved by exploring other test pair correlation functions correponding to disordered packings. We 
demonstrate that there are simpler test functions that lead to the same asymptotic result. More 
importantly, we show that there is a wide class of test functions that lead to precisely the same 
putative exponential improvement and therefore the asymptotic form i/2^^-'^'^^^^---^^ is much more 
general than previously surmised. This class of test functions leads to an optimized average kissing 
number that is controlled by the same asymptotic behavior as the one found in the aforementioned 
paper. 

PACS numbers: 05.20.-y, 61.20.-p 
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I. INTRODUCTION 



A collection of congruent spheres in d- dimensional Euclidean space M.'^ is called a sphere 
packing if no two spheres overlap. Although the practical relevance^ of sphere packings 
in high Euclidean dimensions was appreciated by Shannon in 1948 there has been a 
resurgence of interest in such problems in both the physical and mathematical sciences 



,i3i,Q,fly,0,i8: 
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ll| . Shannon showed that the optimal way of sending digital signals 
over noisy channels corresponds to the densest sphere packing in a high dimensional space [l|. 
These "error-correcting" codes underlie a variety of systems in digital communications and 
storage 



12|. Physicists have investigated sphere packings in high dimensions to gain insight 



into c 



0, y, Q 



assica 
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ground and glassy states of matter as well as phase behavior in lower dimensions 



ITj . Understanding the symmetries and other mathematical properties of the 



densest pac 



dngs in arbitrary dimension is a problem of long-standing interest in discrete 
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geometry |5|: 

The packing density or simply density (/> of a sphere packing is the fraction of space M.'^ 
covered by the spheres. We will call 



sup 0(-P) 



(1) 

the maximal density, where the supremum is taken over all packings in R''. The set of lattice 
packings is a subset of the set of sphere packings in M.'^. A lattice A in M.'^ is a subgroup 
consisting of the integer linear combinations of vectors that constitute a basis for M'^. A 
lattice packing Pl is one in which the centers of nonoverlapping spheres are located at the 
points of A. In a lattice packing, the space M'^ can be geometrically divided into identical 
regions F called fundamental cells, each of which contains the center of just one sphere. 
In the physical sciences, a lattice packing is simply a packings arranged on the sites of a 
Bravais lattice. Non-lattice packings include periodic packings (more than one sphere per 
fundamental cell) as well as disordered packings fisl . 

The sphere packing problem seeks to answer the following question: Among all packings 
of congruent spheres, what is the maximal packing density 0max, i-e., largest fraction of M.'^ 
covered by the spheres, and what are the corresponding arrangements of the spheres {12!. [l^? 
For arbitrary d, the sphere packing problem is notoriously difficult to solve. In the case of 
packings of congruent rf-dimensional spheres, the exact solution is known for the first three 
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space dimensions. For d = 1, the answer is trivial because the spheres tile the space so that 
0max = 1- In two dimensions, the optimal solution is the triangular lattice arrangement 
(also called the hexagonal packing) with 0max = vr/V^. In three dimensions, the Kepler 
conjecture that the face-centered cubic lattice arrangement provides the densest packing 
with 0max = tt/vTs was only recently proved by Hales [51]. For 3 < d < 10, the densest 
known packings of congruent spheres are lattice packings (defined below) . For example, the 
"checkerboard" lattice D^, which is the d-dimensional generalization of the FCC lattice, is 
believed to be optimal in and M^. The and Leech lattices in and M^^, respectively, 
are remarkable dense and symmetric and are most likely the densest packings in these 
dimensions [11]. However, for sufficiently large d, lattice packings are most likely not the 
densest, but it becomes increasingly difficult to find specific dense packing constructions in 



high dimensions 



19\. In high dimensions, the best that one can do theoretically is to devise 



rigorous upper and lower bounds on 



Upper and lower bounds on the maximal density 0max exist in all dimensions [12]. 
Minkowski proved that the maximal density 0max aniong all Bravais lattice packings 
for (i > 2 satisfies the lower bound 



max 



2d-i 

where ({d) = YlT=i Riemann zeta function. One observes that for large values of 

d, the asymptotic behavior of the nonconstructive Minkowski lower bound is controlled by 
2^^^. Since 1905, many extensions and generalizations of ([2]) have been obtained [l2], but 
none of these investigations have been able to improve upon the dominant exponential term 
2"*^. The best currently known rigorous lower bound on ^^g^^ obtained by Ball [0]. He 
found that 

> (3) 

Interestingly, the density of a saturated packing of congruent spheres in for all d satisfies 
the lower bound [22|] 

^ > (4) 
and thus has the same dominant exponential term as the Minkowski lower bound ([2]). A 
saturated packing of congruent spheres of unit diameter and density in M*^ has the property 
that each point in space lies within a unit distance from the center of some sphere. As we 
will discuss below, the lower bound (jlj) is not a stringent bound for a saturated packing and 



4 



therefore is improvable. 



Rogers 
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14| found upper bounds on the maximal density (/)max by an analysis of the 
Voronoi cells. For large d, Rogers' upper bound asymptotically becomes 
batiansky and Levenshtein found an even stronger bound, which in the limit d —>■ oo 
yields 0max < 2~^-^^^^'^^^^°^^^\ Cohn and Elkies obtained and computed linear program- 
ming upper bounds, which provided improvement over Rogers' upper bound for dimensions 
4 through 36. They also conjectured that their approach could be used to prove sharp 
bounds in 8 and 24 dimensions. Indeed, Cohn and Kumar used these techniques to 
prove that the Leech lattice is the unique densest lattice in M^"^. They also proved that no 
sphere packing in M^*^ can exceed the density of the Leech lattice by a factor of more than 
1 + 1.65 X 10~^°, and gave a new proof that the Eg lattice is the unique densest lattice in 

A recent investigation proves that there exists a disordered packing construction in 
with a maximal density that achieves the saturation lower bound (jlj) for any d. This 



construction is referred to as the "ghost" random sequential addition (RSA) packing [24 1 
and it was shown that all of the ri-particle correlation functions for this packing can be 
obtained analytically for all allowable densities and in any dimension. Interestingly, this 
packing is unsaturated (see Fig. 1) and yet it has a maximal density 2"*^, which suggests 
that there exist disordered saturated packings that exceeds the saturation lower bound (jl]) 
or the Minkowski lower bound ([2]). Indeed, another recent study \10\ strongly suggests that 
the standard disordered RSA packing [25] at its maximal saturation density scales as d2~'^ 
for large d, which has the same asymptotic behavior as Ball's lower bound ([3]). Note that 
spheres in both the ghost and standard RSA packings cannot form interparticle contacts, 
which appears to be a crucial attribute to obtain exponential improvement on Minkowski's 
bound [7|, as we discuss below. 

Do there exist disordered packings that can provide the long-sought exponential im- 
provement of Minkowski's lower bound? Torquato and Stillinger [7| employed a conjecture 
concerning the existence of disordered sphere packings and an optimization procedure that 
maximizes the density associated with a "test" pair correlation function to provide the 
putative exponential improvement on Minkowski's 100-year-old bound on 0max (see Sec- 
tion II for details). The asymptotic behavior of the conjectural lower bound is controlled by 
2-{(o.77865...))'^_ Moreover, this lower bound always lies below the density of the densest known 
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FIG. 1: (Color online) A configuration of 468 particles of a ghost RSA packing in at a density 
very near its maximal density of 0.25. This was generated using a Monte Carlo procedure within 
a square fundamental cell under periodic boundary conditions. Note that the packing is clearly 
unsaturated and there are no contacting particles. 

packings for 3 < d < 56, but, for d > 56, it can be larger than the density of the densest 
known arrangements, all of which are ordered. These results counterintuitively suggest that 
the densest packings in sufficiently high dimensions may be disordered rather than periodic, 
implying the existence of disordered classical ground states for some continuous potentials. 
In addition, a decorrelation principle for disordered packings was identified in Ref. [t], which 
states that unconstrained correlations in disordered sphere packings vanish asymptotically 
in high dimensions and that the Qn for any n > 3 can be inferred entirely (up to some 
small error) from a knowledge of the number density p and the pair correlation function 



(72(1") • This decorrelationprinciple 26|, among other things, provides justification for the 
conjecture used in Ref. B. aud is vvidly aud explicitly exhibited by the exactly solvable 
ghost RSA packing process [6(] as well as by computer simulations in high dimensions of the 
maximally random jammed state and the standard RSA packing process [lol ]. 

In this paper, we investigate whether the putative exponential improvement of 



Minkowski's lower bound found in Ref. 



7| can be further improved by exploring other 
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test pair correlation functions. We will show that there are simpler test functions that lead 
to the same asymptotic result. More importantly, we will demonstrate that there is a wide 
class of test functions that lead to the same exponential improvement as in Ref. Q|. 



II. PRELIMINARIES AND OPTIMIZATION PROCEDURE 

A packing of congruent spheres of unit diameter is simply a point process in which any 
pair of points cannot be closer than a unit distance from one another [tJ. A particular 
configuration of a point process in M'^ is described by the "microscopic" density 



oo 

n{r) 

1=1 



J]5(r-x,). (5) 

This distribution can be interpreted in a probabilistic sense which is particularly useful 
for the arguments we will present, even in the limit in which no explicit randomness is 
present, as in the case in which the spheres are arranged on the sites of a (Bravais) lattice. 
We define the n-particle density as the ensemble average 

p„(ri, ...,r„) = / ^ (5(ri - XiJ ... 5(r„ - XjJ \ , (6) 

which is a nonnegative quantity. Henceforth, we will assume that the random process is 
translationally invariant, i.e., statistically homogeneous. It follows that there is no preferred 
origin in the packing and thus the n-particle densities p„(ri2,ri3, . . . , ri„) only depend on 
relative displacements, where rij = rj — ri. In particular, the one-particle density pi(r) = 
(5(r — xi)) = p is a constant called the number (center) density. Note that the packing 
density defined earlier is related to the number density p for spheres of unit diameter via 
the relation 

^ = pv,{l/2). (7) 

where f i(r) = 7r"'/^r"'/r((i/2+l) is the volume of a sphere of radius r. The surface area of such 
a sphere is Si(r) = 27r'^/V'^-Vr(rf/2). If we divide p„ by p", we get the n-particle correlation 
function 5'„(ri2, ri3..., ri„), which clearly is also a nonnegative function. As will become clear 
shortly, the pair correlation function 5'2(ri2) has particular importance to us. If the point 
process is additionally rotationally invariant (i.e., the packing is statistically homogeneous 
and isotropic), the pair correlation function g2{r) depends only on the distance r = |r|. 
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In Ref. {2! , 5f2-invariant processes were examined in order to gain insights about the nature 
of disordered sphere packings. A g2-invariant process is one in which a given nonnegative 
pair correlation 5'2(r) function remains invariant for all r over the range of densities 

< < 0,. (8) 

The terminal density 0^, is the maximum achievable density for the (72-invariant process 
subject to satisfaction of certain necessary conditions on the pair correlation. In particular, 
they considered those "test" g2{r)^s that are distributions on depending only on the 
radial distance r. For any test g2{r) associated with a packing, i.e., g2{r) = for r < 1, they 
maximized the corresponding density 0, i.e., 

max0 (9) 

subject to the following two conditions: 

g2{r) > for all r, (10) 

Sik) = 1 + pi2nf' drr'-' [^^(r) - 1] > for all k. (11) 

Condition ffTTl) states that the structure factor S{k) [trivially related to the Fourier transform 
of g2{r) — 1] must also be nonnegative for all wavenumbers. It is a known necessary condition 



on the existence of a point process 







Recently, Torquato and Stillinger 



J27l|. but it is generally not sufficient [281]. 

conjectured that a disordered sphere packing in M.'^ 
at number density p exists for sufficiently large d if and only if the conditions (fTOl) and 
ffTTj) are satisfied. The maximum achievable density is the terminal density 0*, which then 
implies the lower bound 

0max > 0* (12) 

There is mounting evidence to support this conjecture. First, the aforementioned decorre- 
lation principle states that unconstrained correlations in disordered sphere packings vanish 
asymptotically in high dimensions and that the for any n > 3 can be inferred entirely 
from a knowledge of p and g2- Second, other necessary conditions on g2, such as the Yamada 



condition j29| as well as others |3], appear to only have relevance in very low dimensions. 
Third, one can recover the form of known rigorous bounds [cf. ([2]) and ([3])] for specific 
test 5'2's when the conjecture is invoked. Finally, in these two instances, configurations of 
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disordered sphere packings on the torus have been numerically constructed with such g2 in 



3G, 



3l|. 



low dimensions for densities up to the terminal density 

Interestingly, the optimization problem defined above is the dual of the infinite- 
dimensional linear program devised by Cohn and Elkies 4| to obtain upper bounds on 
the maximal packing density. In particular, let /(r) be a radial function in M'^ such that 

/(r) < for r > 1, 

f{k) > for all k. (13) 
Then the number density p is bounded from above by 

mm — - — . 14) 

2<^/(0) ^ ^ 

The radial function /(r) can be physically interpreted to be a pair potential. The fact 
that its Fourier transform must be nonnegative for all k is a well-known stability condition 
for many-particle systems with pairwise interactions [32-]. We see that whereas the linear 
program specified by flTU]) and flTTl) utilizes information about pair correlations, its dual 
urogram fll3l) and (1141) employs information about pair interactions. It is important to note 
3] that even if there does not exist a sphere packing with g2 satisfying conditions (fTO|l and 
ffTTj) . the terminal density 0* can never exceed the Cohn- Elkies upper bound. Every linear 
program has a dual program and when an optimal solution exists, there is no duality gap 
between the upper bound and lower bound formulations. However, until recently, it was not 
clear how to prove that there was no duality gap for the aforementioned infinite-dimensional 



sphere-packing linear program j3|]. Recently, Cohn and Kumar 33|] have proved that there 
is no duality gap. 

By means of the linear program described above and the aforementioned conjecture con- 
cerning the existence for a certain test function g2, it was found in Ref. j?! that in the limit 
d ^ oo, 

0max > 0* ~ 2-i'^+T^i+2-^2497...dV3+i log, d+log,(3.2761...) ^ (^5) 

where the terms neglected are monotonically decreasing with d. The first term in the series 
provides the putative exponential improvement of Minkowski's lower bound ([2]). In the 
following, we will be interested mainly in the exponential improvement of Minkowski's lower 
bound, and so we simplify the right-hand side of f|T5l) by writing it as 



This is not to be intended as an asymptotic expansion of 0* in the sense of Poincare (the 
ratio of the right-hand side to the left-hand side does not go to unity when d oo), however, 
it is an asymptotic expansion in such sense for logg 0*. 

In what follows, we will show that we can obtain a conjectural lower bound asymptotically 
equal to (1151) with a simpler test function. Then we will demonstrate that the requirement 
of hyperuniformity 27] in Ref. j?] is actually a necessary condition that arises only from 
the optimization procedure. Finally, we will show some examples of how enlarging the 
space of test functions where the optimization is performed does not change the asymptotic 
exponential behavior, although non-exponential improvement is found. 

Although these results do not constitute a proof of lower bounds, they strongly suggest 
that an estimate of the asymptotic behavior of the solutions to the lower-bound linear 
programming problem can be achieved and that physical intuition is gained about the spatial 
structures they describe. 



III. STEP PLUS DELTA FUNCTION REVISITED 

Following Torquato and Stillinger [7|], we choose the following test g2{r): 

g,[r) = e(r - 1) + j^^Sir - 1). (17) 

Here the parameter Z has the interpretation of the average kissing number. The structure 
factor becomes 

'd^ Jd/2-iik ) 

J^d/2-l 

a{k) 2'^(f) + b{k) Z, (18) 



S{k) 



which defines the functions a,b. The terminal density is defined by the linear program (I9l), 
ffTOj) and ffTTj) . Z is then a free parameter to be optimized appropriately. 

Unlike Torquato and Stillinger [3], we do not impose hyperuniformity 27] (requiring the 



structure factor to vanish at A; = 0) to simplify the optimization. Moreover, we are also 
interested in finding the largest average kissing number Z that (for a given d) satisfies the 
constraints. In this latter case, it is that must be chosen appropriately. These are two 
infinite-dimensional, linear programming problems. 

There is a graphical construction that will help us look for such points and that will be 
helpful also in cases where more parameters are to be varied. For any given k the set of 
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allowed points in the {4>,Z) plane [i.e., those for which S{k) > 0] is the half plane above 
(below) the line 1 — a{k)2'^(j) + b{k)Z = for positive (negative) a. Upon changing k by a 
small step to k + A, we repeat the construction and find the intersection of the two half- 
planes. By letting k vary over the positive reals and letting A ^ 0, we find a limiting finite, 
convex region B which gives the allowed values of 0, Z. This region is the set internal to the 
curve obtained by solving the equations 

S{k,^,Z) = 0, ^S{k,(^,Z)=0, (19) 

with respect to 0, Z. This is depicted in Fig|2l It is not difficult to prove that the region B 
is indeed internal to the entire spiral. It will suffice to observe that the distance of a point 
on the spiral from the origin is a monotonically increasing function (for sufficiently large k). 

Now the terminal density 0* is the x-component of the rightmost point in B. Analogously 
the y-component of the topmost point in B gives the terminal kissing number Z^,^,. 

The terminal density is found at the first zero of b{k), which is located at the first zero of 
the Bessel function of order d/2 ~ 1. As customary we call this number jd/2-i,i- The value 
of (0^,, Z^,) is then found by finding the point on the spiral corresponding to k = jd/2,i'- 

J. _ 2-'^ _ g-3rf/2 {jd/2-l,lY^'^ ^20) 



a(jd/2-i,i) r(l + d/2)Jd/2ijd/2-i,i) 

2%. (21) 



a'(jd/2-i,i) _ a'{jd/2-i,i)^d 



b'ijd/2-i,i)a{jd/2-i,i) b'{jd/2-i,i) 
By using the asymptotic formulas, valid for large u 

= u + 1.85576... z/^/^ + O (z/"^/^) , (22) 

J.(>-i,i) = -J'u-iUu^i,i) = 1.11310... u-^/'' + O {u~^'^) , (23) 

we find 

0, ~ 2-i'^+i^i+2.12497...di/3 _ 2-(0.77865...)d_ (24) 

Notice that this is the same case that was treated in [tI but there hyperuniformity was 
imposed and the Minkowski bound was recovered. Here we are not imposing hyperuniformity 
and the resulting terminal structure factor is not hyperuniform. The form of S{k) at the 
terminal point 0*, is given in Figure 2. Notice that the first zero is at A; = jd/2-1,1 — d/2. 
This can be interpreted as the appearance of a structure with length-scale £ ~ l/d in the 
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FIG. 2: (Top panel) For d = 16, the set B of allowed packing densities and kissing numbers. The 
rightmost point is the maximal packing density ^* and its corresponding kissing number Z^. The 
topmost point is the maximal kissing number Z** which corresponds to packing density 0** = 0. 
(Bottom panel) As in top panel, the region B of allowed packing densities and kissing numbers 
for d = 16. For convenience in plotting, the horizontal and vertical axes represent the functions 
e((;/!)) logio(|2'^</'|), and e(Z) log^od-^l), where e{x) = signx, respectively (although in this way the 
small region \2'^<p\ < 1, \Z\ < 1 had to be left out of the graph). This figure shows how the solution 
of the equations S{k,(l),Z) = 0, dS{k,(t),Z)/dk = for varying k form an ever-growing spiral in 
which the allowed region B is completely contained. So this geometrical construction proves that 
every point in B are solutions to the linear programming problem S{k, (p,Z) > 0, (p > 0, Z > 
for every k >0. 12 



CO 




FIG. 3: The structure factor for the terminal density = 0.0046692, = 217.121 in d = 16. 
Notice the zero at k* = jj^i = 11.086... and the rapid asymptoting to the value S{oo) = 1. 

system at large d. However, since a sphere packing corresponding to such an S{k) could not 
be hyper uniform, it cannot be a Bravais lattice. 



Followmg a, we check whether the Yamada condition M on the number variance M 

is satisfied by the pair correlation f[T7l) . As in [7|, we find a violation only for d = 1. 

The terminal kissing number is given by the topmost point in B which is the point k** 
where a{k**) = 0. It can be easily proved that h'{k**) = as well so that (p^^ = and 

1 



2(1172-1)1. 



(25) 



bUd/2,i) 

It is intriguing to notice that the density corresponding to the terminal kissing number is 
zero. 



IV. STEP PLUS DELTA FUNCTION WITH A GAP 

This case was analyzed by 7|] before by imposing hyperuniformity. Here we show that in 
order to find the terminal density, one does not need to impose hyperuniformity from the 
beginning but rather that it arises as a necessary condition form the optimization procedure. 
We will show that the same asymptotic behavior of the terminal density found in the previous 
example is obtained (modulo non-exponential prefactors). 
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We choose the test function 

g,[r) = e(r - (1 + a)) + j^S{r - 1), (26) 

depending on two parameters, Z, a and the density of centers p. Performing the integrals 
gives the corresponding structure factor 

S{k) = 1 - a((l + a)k)2'^{l + aY(P + h{k)Z, (27) 

where the functions a, h were defined in the previous section. Again we look for the rightmost 
point in the set, which is now given by 

(28) 
(29) 



{l + aYa{{l + a)jd/2~i,i) 
{l + a)a'{{l + a)k) 



fc'(jd/2-i,i)a((l + o-)jd/2-i,i) ' 
We now need to maximize the value of 0* over a. Clearly, we can increase a to increase 
0* indefinitely until a((l + o")jd/2-i,i) becomes zero, namely when (1 + cr)jd/2-i,i = jd/2,1, 
which gives a ~ 2/d. The prefactor goes to a constant: (l + cr)"' ~ {1 + 2/dY ~ and 
does not change the asymptotic dependence on d. This would suggest that the density 
can be increased without bound by adjusting the other parameters. This is not the case, 
however, since when we increase a we encounter the first "global" obstacle [by which we 
mean at wavenumbers k far from the first zero of h{k), which was setting the relevant k 
scales up to now] at the value of a when (1 + aY^'^cf)^ = — 1. Notice that a(0) = 
6(0) = 1 and both functions decrease monotonically until their first zeros; here we have 
S{0) = 1 — (1 + aY'^'^4'* + Z^: = and any further increase of a would make 5(0) < 0. Thus, 
hyperuniformity has arisen as an optimality condition. Of course one should make sure that 
there is not a disconnected region in the parameter space {a, (p, Z) with better terminal 
density 0* but where hyperuniformity does not hold. We have searched the parameter space 
by discretizing the relevant range of k and solving, using Mathematica, the linear program 
&, (ITOl) and ( ITTl) . We have not been able to find another allowed region of the parameters 
disconnected from the previous one. 

Hence we assume that the global terminal value 0* is indeed obtained by imposing hype- 
runiformity and maximizing with respect to the remaining parameters (the two operations 
can be performed in any order). 
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We notice that now we have reduced the problem to the case that has been analyzed in 
Q]. We will not repeat that analysis here but refer the reader to that paper. It is important 
to observe that in 7j] the resultant asymptotic scaling law for the terminal fraction 0* 
coincides with the one presented in the previous section 0^, ~ 2~(*'-'^''^^^ - )'^. Although the 
non-exponential terms are different from those in the previous section, it is remarkable that 
the same exponential scaling law arises for two different cases. This strongly suggests that a 
large class of test functions can possess this asymptotic behavior. With this in mind, we go 
on to analyze the next case win which the test pair correlation function consists of a hard 
core with two delta functions and a gap. 

V. STEP PLUS TWO DELTA FUNCTIONS WITH A GAP 

In this Section we find the solution of the optimization problem (Q, (|TOl) and ( fTTl) for the 
family of pair correlation functions g2{r) composed of unit step function plus a gap and two 
delta functions, one at contact and the other at the end of the gap: 

g,{r) = e{r - (1 + a)) + -^6{r - 1) + 6{r - (1 + a)). (30) 

s{l)p s{l + a)p 

This family depends on 3 parameters, a, Zi, Z2 and we need to optimize them in order to 
find the optimal terminal density 0*. The structure factor is 

- W/2+l)(l + .)^2-/^^||^^±^ (31) 

= I + Z2 c{k) + Zi b{k) - (1 + aY2'^(j) a{k), (32) 

where the last line defines the functions a, b, c. Notice that a(0) = 6(0) = c(0) = 1 and 
|a(A;)|, |6(/i;)|, \c{k)\ < 1 follow from the properties of the Bessel functions. It is also con- 
venient to reabsorb the factor (1 + cr)°'2'^ in the definition of 0, i.e. (1 + cr)°'2°'0 — > 0. We 
will restore the proper units at the end of the calculation. The solution of this optimization 
problem for arbitrary d is a formidable task. However, guided by the results of the previous 
section, we assume we can find an improvement on the previous bound even after imposing 
hyperuniformity. 

Therefore, we fix the value of Z2 = (f) — Zi — 1 and are left with the other two parameters 
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to optimize. Inserting this value of in (1321) . we find the reduced optimization problem 



S{k) = (1 - c{k)) - [a{k) - c{k))<\) + (6(A;) - c{k))Zx > 0. (33) 
By using the fact that c{k) < 1 we might as well study the optimization problem 

S^'\k, a, 0, Zi) = = 1 - a{k)(^ + > 0, (34) 

- (35, 

^ M. ,3a) 

Formally, this problem is analogous to the previous case with one delta function with gap 
and can be studied in the very same fashion. The process of having solved for Z2 and 
changed the functions a, 6 to a, /5 can be thought of as a renormalization process that allows 
to integrate out one delta function to reduce the problem to a simpler one. 

The mathematical problem of finding the terminal fraction is formally identical to that of 
the previous section, although the constitutive functions a, (3 are more complicated. How- 
ever, as long as a numerical analysis is concerned this does not present further difficulties. 

We proceed in the following way: for a fixed cr we find the rightmost point of allowed 
region, 0*(o"), ^1^^,(0"), by finding the first zero of P{k), call it k*, 

, . . 1 



cr 



a{k* 



(37) 



We then maximize the value of 0*(cr) with respect to variations of a. Generically, in- 
creasing a increases the value of 0* until a positivity condition is violated (for small k). It 
turns out that the first condition to be violated is 5'*^^''(0) > 0. So in practice we find the 
terminal value of a by solving the equation 

5«(0,fx,</)*(a),Zi,,(a)) =0, (39) 

with respect to a. Notice that this is now a "strong" hyperuniformity requirement, since 
S^^\k) ~ fc^ near the origin implies S{k) ~ k^ near the origin, since 1 — c{k) ~ k'^. We 
are tempted to conjecture that this is a universal feature: adding more delta-functions to g2 
and solving the linear programming problem, we obtain structure factors S{k) that become 
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increasingly flatter at the origin. Hence, at least in this respect, the structure factor looks 
increasingly similar to that of a lattice. 

As can be seen from Table 1 and Figure H] (here the proper normalization for has been 
restored) the improvement on the previous bound is relevant but the asymptotic exponent is 
the same. Analytically, it is not difficult to obtain the rate of exponential decay (dictated 
mainly by the Stirling expansion of the gamma functions and the scaling of the first zero of 
(3 with d for large d), which turns out to be the same as the previous cases, namely. 

It is plausible, therefore, that the incorporation of any finite number of delta functions in a 
test g2 will not improve the exponent in fHOj) . This exponent fits the numerical data very 



well. A best fit of the data in Table HTl using the functions d, d^^^ ,\og2d, appearing in the 
analysis in the previous section and invoking the existence conjecture of Ref. [3] yields the 
putative lower bound 

The first term is fixed by our analysis, the d^^^ is term is consistent with the analytic value 
2.1247 in Eq. ( JT5l) . The sub-leading term logj d in this expression is very difficult to obtain 
analytically and we have not succeeded in this task. However, it is clear that there there is 
an improvement from the value 1 = 0.1666 . . . appearing in (|T5l) . The improvement is also 
evident from the numbers in Table [B 

It is worth noting that for large d the optimum gap a ~ ^''^J'" (from a best fit analysis). 
This scaling with d is slightly different from that found in the previous section and in 71] 
(there a ~ 1.81/d). Again notice that the scaling of a with d, a oc 1/d is necessary in order 
not to introduce an exponential suppression of density. In fact for large d, (1 + c/dY e'^ 
multiplies the density in all the formulas (and hence it reduces the terminal value by e~'^) . 
A larger gap, say 0{d^^^^^^), would suppress the density by an exponentially large amount 



Table [T] compares the final results of our analysis for the conjectured lower bound on the 
maximal density to the previous lower bound, the best known packings, and the optimal 
upper bound in {4] for selected dimensions up to (i = 150. As in the previous cases, the 



Yamada condition 



2^ is violated only for d = 1. This supports the conclusion reached in 
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Ref. j7| that the Yamada condition appears to only have practical relevance in very low 
dimensions. 



TABLE I: Estimates of the maximal densities for selected dimensions up to d = 150. (/)f,,fc is the 
densest known packing, (f)cE is the upper bound of Cohn and Elkies, (pi^^, is the terminal density 
for a single delta function and 02,* for two delta functions. 



d 








V^*,2 


3 


0.74049 


77982 


U.D (ODD 




63306 


4 


61685 


0.64774 


U.4Z0/D 




47885 


5 


n 46527 


52506 


u.ouoyi 




35437 


5 


37295 


0.41776 


U.ZioDU 




24966 


7 


29530 


32757 


U.i4 ( io 




0.17991 


« 


n 25367 


95367 


U.Uyyoo 




19467 


12 


n 04945 


08384 


U.Uiyio 




095721 


1 ^ 

-LtJ 


01 685 


03433 


U.UUOiD 




00792 


19 


004121 


009885 


000845 




001233 


24 


0.00193 


0.00193 


8.24 X 10" 


-5 


0.000125 


31 


1.18 X 10"^ 


1.93 X 10""^ 


2.91 X 10" 


-6 


4.57 X 10"^ 


36 


6.14 X lO"'^ 


3.59 X 10"^ 


2.57 X 10" 


-7 


4.13 X lO"'^ 


56 


2.33 X 10^1^ 




1.25 X 10" 


11 


2.13 X 10"^^ 


60 


2.97 X 10"^^ 




1.67 X 10" 


12 


2.87 X 10"^2 


64 


1.33 X 10"^^ 




2.22 X 10" 


13 


3.83 X 10"i3 


80 


1.12 X 10"^^ 




6.52 X 10" 


17 


1.15 X 10-^^ 


100 






2.28 X 10" 


21 


4.11 X 10"2i 


150 


8.44 X 10"^^ 




1.27 X 10" 


32 


2.30 X 10"^^ 



VI. CONCLUSIONS AND OPEN QUESTIONS 

The problem of finding the asymptotic behavior of the maximal density (/)max of sphere 
packings in high dimensions is one of the most fascinating and challenging problems in 
geometry. In this paper, we have shown how, using linear programming bounds and a 
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Space Dimension, d 



FIG. 4: (Color online) Comparison of different results for the maximal density 0max versus di- 
d. From bottom to top: Torquato-Stillinger result \h\ (1-delta function with gap), one 



mension 



of the results of this paper (2-delta functions with a gap), densest known packings 
Cohn-Elkies upper bound [4]. 



]J], and the 



TABLE II: Terminal density 0* for two delta functions and a gap, corresponding optimal gap cj, 
and optimal average kissing number Z\ ^, for large d. 



d 







0* 


200 


0.013508 


1.57 X 10^^ 


1.06 X 10^ 


-43 


250 


0.010895 


7.15 X 10^1 


4.18 X 10^ 


-55 


300 


0.009132 


2.94 X 10^5 


1.49 X 10^ 


-66 


350 


0.007862 


1.12 X 10^9 


4.96 X 10^ 


-78 


400 


0.006903 


2.93 X 10^5 


1.56 X 10" 


-89 


450 


0.006154 


1.38 X 10^^ 


4.73 X 10" 


101 


500 


0.005553 


4.67 X 10^^ 


1.40 X 10^ 


112 



conjecture concerning the existence of disordered sphere packings based on pair-correlation 
information, the asymptotic conjectural lower bound 7| 



- > o-(0.77865...)d 



(42) 
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which provides the putative exponential improvement on Minkowski's century-old lower 
bound ([2]), is actually much more general than one could have initially surmised. Precisely 
the same exponential improvement arises for a simpler pair-correlation function than the 
one employed in [7] and survives also to a considerable enlargement of the family of test 
functions g2- This family of functions includes two delta functions with a gap (which we 
have shown improves upon the prefactor multiplying 2~(''-'^''^^^ - )'^ given in Ref. [7|]) and, we 
argue, any finite number of delta functions. If this is true, as we believe, it signifies that the 
decorrelation principle alone has a huge predictive power, since an exponential improvement 
of Minkowski's bound has proved to be an extremely difficult problem. 

One outstanding open question is certainly in which sense this is to be interpreted as 
an asymptotic bound. Based on our present, limited knowledge of optimal sphere packings, 
we foresee diverse scenarios. In one case, for sufficiently large d, the importance of higher- 
order correlations is to be neglected altogether and the bound becomes exact by virtue of 
the decorrelation principle. This would mean that the asymptotic Kabatiansky-Levenshtein 
upper bound is far from optimal: a provocative possibility. In a second scenario, it could 
be that "special dimensions" continue to exist for which the negligence of higher-order 
correlations is impossible. In this case, the lower bound obtained by our methods would not 
apply to these special dimensions but will continue to apply to the other dimensions. On the 
other hand, if the frequency of appearance of these dimensions over the integers is decreasing 
then the decorrelation principle is safe. A third but more pessimistic possibility is that these 
dimensions are actually becoming more and more frequent, and our conjectural bound would 
apply only to the subset of dimensions remaining. However, there is absolutely no evidence 
at present for either the second or third scenario. Our best guess at the moment is that the 
optimal packings in very high dimensions will possess no symmetry at all and therefore are 
truly disordered. If so, then the decorrelation principle dictates that pair correlations alone 
completely characterize the packing in high d, implying that the form of the asymptotic 
bound (H2I) is exact! 

The fact that pair correlations can completely specify an optimal packing may seem to 
be counterintuitive at first glance, but we can now identify even low dimensions where this 
phenomenon occurs. Specifically, whenever the linear programming bounds are exact (i.e., 
achieve some packing), pair correlation information is sufficient to determine the optimal 



packing! This outcome, in all likelihood, occurs in M , M and M lUl]. This implies 
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that whenever hnear programming; bounds are not sharp in low dimensions (albeit without 
a duality gap for any d [33.]), information about high-order correlations are required to get 
optimal solutions. 

Another interesting question arises because our procedure, like Minkowski's, is noncon- 
structive. Specifically, it is an open question whether there exist packing constructions that 
realize our test (72's. For future investigations, it would be fruitful to determine whether 
there are periodic or truly disordered packings that have pair correlation functions that 
approximate well the ones studied in this paper. If these packings could be identified, one 
should attempt to ascertain whether the higher-order correlations diminish in importance as 
(i — > 00 in accordance with the decorrelation principle. If such packings exist (or better, if a 
(i-dependent family of them does), they would enable one to place on firm, solid ground the 
putative exponential improvement on Minkowski's bound. We are currently investigating 
these questions. 
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